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Estimating Channel Impulse Response and Equalizer Coefficients in 
UWB Conmiunication Systeriis 

Cross-Reference to Related Applications 

[01] A claim of priority is made to U.S. Provisional Patent Application 
Serial No. 60/45 1 ,528, filed March 3, 2003, entitled "Channel Estimation 
and Equalizer Coefficient Estimation for UWB Radio Communication 
Systems." 

Field of the Invention 

[02] The present invention relates generally to radio communication 
systems, and more particularly to channel and equalization estimation in an 
ultra wide bandwidth conmiunications systems. 

Background of the Invention 

[03] With the release of the "First Report and Order," Feb, 14th, 2002, 
by the Federal Communications Commission (FCC), interest in ultra wide 
bandwidth (UWB) communication systems has increased. The IEEE 802.15 
standards organization, which is responsible for Personal Area Networks, 
has established a task group, TG3a, to standardize a high-data-rate physical 
layer based on UWB . 

[04] Ultra wide bandwidth (UWB) communication systems transmit 
and receive extremely short electro-magnetic energy impulses, therefore the 
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terms "UWB" and "impulse radio" are used synonymously herein. Because 
the bandwidth of the pulses is much wider than the bandwidth of the payload 
signal, UWB is a form of spread-spectrum communication. Each pulse can 
cover anywhere from 500 MHz to several GHz of the radio spectmm. 

[05] Information is transmitted by modulating the frequency, timing, 
amplitude, polarity, or other aspect of the pulses. UWB systems can provide 
very high data rates for short-range wireless communications. In part, UWB 
systems are designed to distribute information in home, office, school, and 
industrial environments using high-speed links between computer devices. 

[06] However, it is a challenge to realize UWB modulation schemes 
that increase data throughput while minimizing errors in dense multi-path 
environments. Propagation measurements and channel modeling studies 
indicate that many multi-path components (MFCs) are expected for an UWB 
radio link. In a dense multi-path environment, the number of MFCs 
increases linearly with the bandwidth. 

[07] For example, a system with a 10 GHz bandwidth, operating in an 
environment with a maximum excess delay spread of 100 ns requires a 
thousand fingers in a rake receiver. Even in a sparse multi-path environment, 
like the channels specified by the IEEE 802. 15.3a standard for channel 
models, up to 80 fingers are required to collect 80% of the available energy. 

[08] In order to collect most of the available energy from MFCs, while 
at the same time reducing the number, of fingers in the rake receiver, it is 
desired to estimate the channel impulse response so that the correct reference 



2 



MH-5144 
Molischetal. 

waveforms for convolutions with the received signal at each rake finger can 
be determined. 

[09] The requirements for channel estimation in the IEEE 802.15.3a 
standard for data format are stringent. Data are to be transmitted in blocks of 
8,000 bits. It is also assumed that the channel can change from block to 
block. Each block of data lasts only 0.1 ms. The typical speed of movement 
in an indoor environment is about 1 meter per second, so for a 10 GHz upper 
frequency limit, the maximum Doppler frequency is 30 Hz. Thus, the 
channel stays stationary for only 30 ms. 

[010] The impulse responses can last up to 200 ns, as can be seen in the 
Final Report and Order. Because sampling has to be done at a rate of about 
10 G samples per second, this means that 2000 samples of the channel 
impulse response have to be evaluated. If data are transmitted in a constant 
stream, then is should be possible to exploit the correlation between channel 
realizations. 

[Oil ] In the prior art, the channel impulse response and equalizer 
coefficient are usually estimated from a single pseudo-noise (PN) training 
sequence. However, other users and out-of-band interferers might be 
transmitting while the single training sequence is transmitted. Therefore, 
there is a need to suppress co-channel interferers. This can be done by 
spreading the training sequence. 

[012] A brute-force approach samples and A/D converts the measured 
impulse response at a speed of 10 G samples per second. In principle, the 
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channel sounding requires only two microseconds, i.e., the length of the 
impulse response times a factor of 10 for the interference suppression, which 
is 200 symbol durations. This is much shorter than the duration of the data 
block, and thus is not a significant overhead. However, A/D converters that 
can process 10,000,000,000 samples per second are prohibitively expensive. 

[013] During the estimation period, while the single training sequence is 
transmitted, it is necessary to sample the received signal at a chip rate to find 
the delays and amplitudes of the received multi-path components. Because 
the number of MFCs is not yet known at this time, each possible resolvable 
delay, i.e., each chip period, must be estimated. For a 10 GHz system, that 
again means sampling and A/D converting at 10 G samples/s, which is 
prohibitively expensive. 

[014] Therefore, there is a need for method and system that can estimate 
the channel impulse response and equalizer coefficient at a lower cost 
without degrading performance. 

Summary of the Inveiition 

[015] It is an object of the invention to provide a low cost channel 
estiniator that can estimate the channel impulse response and equalizer 
coefficients in an ultra wide bandwidth (UWB) communications system. 

[016] According to the present invention, multiple training sequences are 
transmitted. Each training sequence includes, for example, fifteen symbols, 
and the training sequence is repeated, for example, fourteen times to 
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suppress interference. There is a guard interval between two consecutive 
training sequences to avoid inter-sequence interference. Each symbol is 
'spread' over the spectrum by a predetermined sequence of chips at a chip 
rate. 

[017] After front-end processing, the UWB signal of the training 
sequences is passed to a rake-like receiver structure with a small number of 
parallel correlators (fingers), for example, three, for channel impulse 
response estimation. 

[018] The fingers sample each received training sequence, in parallel, 
with different time offset to produce samples. The samples substantially 
span the entire time interval of the impulse response. If the training sequence 
is repeated, then the samples have a fine resolution, and the channel impulse 
response can be estimated from the transmitted signal and the samples. 

[019] After the channel has been estimated, the weights for the rake 
receiver used for the data reception are determined. Finally, equ^izer 
coefficients are estimated by minimizing the error between equalized outputs 
and the equalizer training sequence. 

Brief Description of the Drawings 

[020] Figure 1 is a detailed block diagram of a receiver with channel 
impulse response and equalizer coefficients estimation according to the 
invention; 
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[021] Figure 2 is a timing diagram of multiple training sequences 
according to the invention; 

[022] Figure 3 is a flow diagrani of a method for the estimating 
according to the invention; 

[023] Figure 4 is a graph comparing results obtained by the estimating 
according to the invention; 

[024] Figure 5A is a graph showing bit error rates; 

[025] Figure 5B is a graph showing normalized mean-square errors; 

[026] Figure 6 is a block diagram of a UWB system according to the 
invention; and 

[027] Figure 7 is a block diagram of a UWB receiver according to the 
invention. 

Detailed Description of the Preferred Embodiment 
System Structure and Operation 

[028] Figure 6 shows a UWB system according to the invention. The 
system includes a transmitter 601 and a rake-like receiver 100 connected via 
a wireless channel 603. Both the transmitter and the receiver include a clock 
(CLCK) 604. The clocks are synchronized 605 to each other. The details of 
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how this is achieved can vary, see, e.g., Meyr et aL, Synchronization, 
Channel Estimationy and Signal Processings Volume 2, Digital 
Communication Receivers, Wiley, 1997. The transmitter sends m training 
sequences 200 as described in detail below. The pulses of the training 
sequences are modulated at a chip rate. 

[029] As shown in Figure 7, each of the m training sequences 200 is 
passed^ in parallel, through multiple correlators 701 (rake fingers) of the 
receiver 100. Parallel "sample and hold" circuits 702 sample each received 
training circuits at multiple delay points with fixed offsets to obtain at least 
one (^==1) sample for each training sequence. It should be noted that each of 
the n fingers can also obtain serially more samples (k>\). Thus, the total 
number of samples can be expressed as (m x n x A:). The idea is that these 
samples span a tinie interval corresponding to the impulse response at a 
resolution that is substantially equal to the chip rate. 

[030] However, the actual sampling rate according to the invention is at a 
substantially lower rate then the prior art sampling at the chip rate (on the 
order of 10 GHz) because the siampling extends over a longer period of time 
of the multiple training sequences. In fact, the sampling rate can be at the 
symbol rate (1 lOMHz in an 802.15.3a system). It should be understood that 
the sampling can also be at other low sampling rates. Thus, low cost 
components can be used for the receiver according to the invention, and the 
primary object of the invention is achieved. 
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[031] After the sampling, the analog signal is A/D converted 703 to a 
digital signal. The (mxnxk) digital signals from the fingers are further 
processed 705 to determine the channel impulse response. 

[032] In a second, optional, step, the estimate of the channel impulse 
response is then used to determine the weights of the rake receiver. For the 
basic prior art method see Molisch, Wideband Wireless Digital 
Communications, Prentice-Hall, 200 1 . 

[033] The ralce coefficients are then used to weigh the output of the 
correlators 701 so that most of the available energy of the signal can be 
acquired. Thus, the invention estimates an "effective channel" from the 
antenna of the transmitter to the output of the summation device. Optionally, 
additional training sequences can be sampled at the symbol rate to provide 
an estimate of the channel impulse response of the effective channel, which, 
in turn, is the basis for the determination of the coefficients of an equalizer. 
The low rate samples are thus used to determine a full-resolution impulse 
response of the channel, and equalizer coefficients for the receiver. All of 
this is now described in greater detail. 

[034] Signaling Format 

[035] In the transmitter, symbols are generated at a symbol rate l/T^. 
Transmitted pulses p corresponding the symbols are in the form of a 

Gaussian function p{t)= K^e , whose Fourier transform is 
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p{q)) = K^e ^ . In this example, the variance cP' is determined by 

\\P{0)fdQ) 

^ — , where /is the frequency. From p{t\ the spreading waveform 

\\P(p)fd(0 

can be obtained by w{t) = YjiS^p{t-kT^), where lITc is the chip rate, which is 

Jt=0 

substantially faster than the symbol rate. An example spreading sequence is 
[si^] = {-1 -1 -1 +1 +1 -1 +1+1}. If the symbol sequence to be transmitted is 

{bk}, then the transmitted signal can be expressed as s(t) = ^b^y^{t - kT,) > 
[036] Detailed Receiver Structure 

[037] Figure 1 shows the rake-like receiver 100 according to the 
invention in greater detail. The receiver includes a conventional front-end 
101 and three fingers 110 to acquire most of the available energy of a 
received signal 102, Each finger 1 10 includes a multiplier and low-pass filter 
1 1 1, a programmable training pulse generator 112. These comprise the 
correlators 701 and sample and hold circuits 702 of Figure 7. 

[038] Adjustable weight blocks 113 operate on digital signals 1 14. For 
clarity the details of the sampling and A/D circuits of Figure are not shown 
here. The digital signal is functionally related to a product of the output of 
the front end and the output of the programmable pulse generator. 
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[039] The receiver also includes a timing controller 140 for the 
programmable training pulse generators 112. A channel estimator 150 
provides input for the timing controller and an estimate 151 of the channel 
impulse response, which is used by the timing controller 140. A weight 
controller 155 adjusts the weights of the weight blocks 1 13 of each finger 
according to the estimated channel impulse response. An equalizer (EQ) 
training sequence generator 160 generates an equalizer training sequence 
161 for an equalizer estimator 170 that produces coefficients 171 for a three- 
tap equalizer 130. The receiver produces output 131 is also fed back to the 
equalizer estimator. 

[040] After summing 120, a combined output 121 can be expressed as 
. , where r{t)= js{t + T)^*{T)dT, Xk is the 

impulse response of the channel at rioA, where A is a minimum timing offset 
bet^veen two fingers. 

[041] To acquire as most of the available energy, we select {ni,n2, its} 
such that [y^^ ,Yr,, > Yn, } correspond to the three taps with largest energy 
magnitudes. 

[042] The three-tap linear equalizer 130 is used to mitigate residual inter- 
symbol-interference (ISI). If the coefficients of the equalizer are {c_pCo,Ci}, 

then the equalizer output 131 is 6„ = Y . 

\ 
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[043] To optimize performance, the equalizer coefficients 171 are 
selected to minimize a mean square-error of its output, that is 

MSE^E\b-b^ . ■ ' 

[044] Channel Estimation and Equalizer Coefficient Estimation 
[045] Training Sequences 

[046] Figure 2 shows the structure of the m pseudo-noise (PN) training 
sequences 200 according to the invention. Each of the multiple training 
sequences uses fifteen symbols 201-215- Each symbol is 75 jisec in4^ 
fij^w^d^ If the 

sequence is repeated fourteen times 21 1-224, then the channel coefficients 
can be estimated in 17.5 |ULsec, followed by 2.5 |LLsec for equalizer 
coefficients estimation, for a total of 20 (Xsec. The sequences can be the 
same or different. 

[047] Channel Impulse Response Estimation 

[048] Each training sequences are sampled, in parallel, at a symbol rate 
(l/Ts). During each symbol period, we observe at least n samples, e.g., three 
if there are three fingers. Because we need to estimate the entire impulse 
response at a resolution of the chip rate, we need forty samples during the 
symbol interval. Therefore, we repeat the training sequence fourteen (40/3) 
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times. The timing offset of each finger ean be adjusted by the controller 140 
so that the output of each finger is as indicated in Table A below. 



Table A 



Sequence 


Output of the 1st 


Output of the 


Output of the 


No. 


flnger 


2"'' finger 


3*^** finger 


1 


r(40nA) 


r(40nA + 14A) 


r(40nA + 28A) 


2 


r(40nA + A) 


r(40nA + 15A) 


r(40nA + 29A) 


3 


r(40nA + 2A) 


r(40nA + 16A) 


r(40nA + 30A) 


4 


r(40nA + 3A) 


r(40«A + 17A) 


r(40nA + 31A) 


5 


r(40nA + 4A) 


r(40nA + 18A) 


r(40nA + 32A) 


6 


r(40nA + 5A) 


r(40/iA + 19A) 


r(40nA + 33A) 


7 


r(40nA + 6A)- 


/■(40nA + 20A) 


r(40nA + 34A) 


8 


r(40nA + 7A) 


r(40nA + 21A). 


r(40nA + 35A) 


9 


r(40nA+8A) 


r(40/iA + 22A) 


r(40nA + 36A) 


10 


r(40nA + 9A) 


r(40nA + 23A) 


r(40«A + 37A) 


11 , 


r(40nA + 10A) 


r(40nA + 24A) 


r(40nA + 38A) 


12 


r(40nA + llA) 


r(40nA + 25A) 


r(40nA + 39A) 


13 


r(40nA + 12A) 


r(40nA + 26A) 




14 


r(40nA + 13A) 


r(40«A + 27A) 





[049] If the training sequence 6f/ s for A: = 0, 1 , . . . , 14, then the training 

14 

sequence can be expressed as st{t) = ^btf^w{t-kT^), and the channel output is 

Jt=0. 
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xt{t) = h{t)* st{t)-\' n{t), where h(t) is the channel's impulse response and n{t) is 
additive white Gaussian noise. 

[050] The A-sampled output of each correlator is 

r{n/s) = jjcr(nA + r)st*{T)dT . If the training sequence st(t) is white, then the 

channel parameters can be estimated by h{nA) -^7^- However, st(t) is not 



perfectly white. Therefore, we modify the above estimation by the following 
discrete Fourier transform (DFT) steps: 

1) Find the DFT of h(n^) , //, = DFT{h(nA)}. 

2) Find the DFT of stit), St, = DFT{st{nA)} . 

3) Modify the DFT of M«A) by 

4) Find the modified channel estimation by h{n&) = / - DFT{H, } , 
[051] Improved Channel Estimation 

[052] Because the channel impulse response is correlated for the sub- 
channels, the estimation of the channel impulse response can be improved, 
as shown in Figure 3. For the first few, e.g., two, training sequences, set the 
timing offsets of the rake fingers of different sub-channels for rough channel 
estimation 310. This locates the portions of the impulse response that 
contain most of the available energy. Identify 320 the largest (six) taps 
according to the rough channel estimation. That is, set the offset timing so as 
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to only sample near the rough estimates, instead of over the entire span of 
the impulse response. Then, estimate 340 the impulse response more 
accurately using another four training periods or 5000 nsec, and identify 350 
the largest three taps from the six taps for the rake receiver and channel 
equalization. 

[053] As shown in Figure 8, the rough estimates 801 for an impulse 
response 800, are used as a basis for the accurate estimates 802. Obviously, 
the rough to accurate refinement can be done in any number of steps. 

[054] It should be noted that more then one sample can be acquired 
concurrently by each rake finger for each training sequence. For example, if 
k (two or more) samples are taken at different time offsets, for each 
sequence, then the number of required sequences can reduced. 

[055] This simplified estimator only needs six training sequences instead 
of fourteen, which only requires 7500 nsec. Together with another 2500 nsec 
for equalizer coefficient estimation, the overall simplified estimation takes is 
only 10 [xsec, which is half the time required to estimate the channel from 
fourteen sequences. 

[056] It should also be noted, that the channel is re-estimated for each 
transmitted block. If the channel is slow- varying, a previous estimate can be 
used as a starting point for a next estimate to further reduce the required 
time. 
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[057] Equalization Coefficient Estimation 

[058] After the impulse response is estimated, the rake receiver can be 
used to acquire most of the available energy of the received signal, and to 
mitigate the residual ISI. Therefore, after processing the training sequences 
for channel estimation, another training sequence is sent for equalizer 
coefficient estimation. 



[059] If the output of the rake receiver is x[n] when the training sequence 
is ckk., then the equalizer coefficients can be estimated by minimizing 
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^c,x[n-k]^ct^ 



'x[n + l]] 

x[n] 
x[n-l]j 



Consequently, 
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1 ^ f 1 

TtZ^*'^^^ , where 



[060] Figure 4 shows results for the channel estimation according to the 
invention. When the SNR is below 0 dB, the normalized MSE (NMSE) of 
channel estimation decreases with an increase of the SNR. However, when 
SNR is above 0 dB, the NMSE has a floor, which is because training signal 
is not exciting enough. 



[061] The modified approach can significantly improve estimation 
performance. It can improve NMSE of overall channel estimation (whole) 
by 2 dB and NMSE of the three greatest taps (3-GT) by 2.6 dB . 
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[062] For the three-tap equalizer output according to the invention, 
Figure 5A shows the bit error rate (BER), and Figure 5 A shows the 
normalized NMSEe. For a 1% BER, the required SNR is about -2.5 dB. 

[063] Although the invention has been described by way of examples of 
preferred embodiments, it is to be understood that various other adaptations 
and modifications can be made within the spirit and scope of the invention. 
Therefore, it is the object of the appended claims to cover all such variations 
and modifications as come within the true spirit and scope of the invention. 
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